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Abstract 

We compare two theoretical approaches to the data analysis of the Cassini relativity 
experiment based on the Doppler tracking and the time delay technique that were 
published correspondingly by Kopeikin et al in Phys. Lett. A 367, 276 (2007) and 
by Bertotti et al in Class. Quant. Grav. 25, 045013 (2008). Bertotti et al believed 
that they found a discrepancy with our paper and claimed that our analysis was 
erroneous. The present paper elucidates, however, that the discrepancy is illusory 
and does not exist. The two techniques give the same result making it evident that 
the numerical value of the PPN parameter 7 measured in the Cassini experiment is 
indeed affected by the orbital motion of the Sun around the barycenter of the solar 
system. 
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In 2002 a measurement of the effect of solar gravity upon the phase of coherent 
microwave beams passing near the Sun has been carried out with the Cassini 
mission, allowing a very accurate measurement of the PPN parameter 7 [TITS] . 
The data have been analyzed with NASA's Orbit Determination Program 
(ODP). Relativistic ranging time delay, incorporated to the NASA ODP code, 
was originally calculated by Moyer [3] under assumption that the gravitating 
body that deflects light, does not move. Regarding the Sun, it means that 
the ODP code derives the ranging delay in the heliocentric frame. Let us 
introduce the heliocentric coordinates X a = (X°,X l ) = (cT,X), and use 
notation x a = = (ct,x) for the barycentric coordinates of the solar 

system, which origin is at the center of mass of the solar system. The Sun moves 
with respect to the barycentric frame with velocity v & = dx®/dt amounting 
to ~ 15 m/s due to the cumulative gravitational attraction of Jupiter, Saturn, 
and other planets [I]. We have discovered [5] that though this velocity looks 
small, it affects the measured value of the PPN parameter 7 and can not be 
neglected in the data analysis of such high-precision relativity experiment as 
Cassini [112] . 

A legitimate question arises whether the ODP code accounts for the solar 
motion or not. We analyzed this question in [5] by doing calculations of the 
Doppler shift caused by the gravitational field of the moving Sun. We came to 
the conclusion that the ODP code does take into account the motion of the Sun 
and that this motion affects the Cassini ranging data. The original papers on 
Cassini experiment [Tj2j did not analyze the impact of the solar motion on the 
results of the experiment. Hence, Bertotti et al [6] have decided to repeat our 
analysis [5] by making direct use of a different approach based on the light-time 
equation that is the Shapiro ranging time delay. Unfortunately, insufficiently 
elaborated comparison of the two different mathematical techniques did not 
allow Bertotti et al [6j to reproduce our results. The goal of the present letter 
is to show that the result of Bertotti et al paper (6] exactly coincides with that 
obtained earlier in our paper [5]. 

The ranging time delay in the heliocentric coordinates with the Sun located at 
the origin of this frame is well-known [7]. After making use of the heliocentric 
coordinates it reads 
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where 7 is a parameter of the PPN formalism [7j, X2 and X\ are the helio- 
centric coordinates of observer on the Earth and emitter (Cassini spacecraft) 
respectively, distance of the emitter from the Sun is R2 = \X 2 \, distance of 
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the observer from the Sun is Ri = X\, and R± 2 = \X 2 — Xi\ is the null-cone 
heliocentric distance between the emitter and observer. This equation coin- 
cides exactly (after reconciling our and Moyer's notations for distances) with 
the ODP time-delay equation (8-38) given in section 8 of the ODP manual on 
page 8-19 [3]. 

Moyer [3] had transformed the argument of the logarithm in the heliocentric 
ranging delay (J2J) to the barycentric frame by making use of substitutions 



X 2 => r 2 = x 2 - x Q (t 2 ) 



X 1 ^r 2 = x 1 - sb (*i 



(3) 



where x 2 = x(t 2 ), x% = x(ti) are the barycentric coordinates of the observer 
and the emitter taken at time of observation, t 2 , and emission, t\, respectively. 
The ODP manual [3] does not provide any evidence that these substitutions in 
the ranging time delay (J2J) are consistent with relativity and do not violate the 
Lorentz symmetry of the Cassini experiment. Nonetheless, equations (j2j), ([3]) 
are legitimate transformations from the heliocentric to the barycentric frame 
in the sense that they take into account velocity of the Sun in the argument 
of the ranging time delay in the linearized approximation. This is because the 
solar barycentric coordinate x & (ti) at the time of emission t\ is not the same 
as the solar coordinate x Q (t 2 ) at the time of observation t 2 



x Q {t 2 ) = aj (t!) + w (ti)(t 2 - h) + O - t^ 2 ) 
due to the non-zero velocity v Q (ti) of the Sun. 



(4) 



However, spatial transformations are not sufficient in order to get all 
velocity-dependent terms of the first order in the ranging time delay. The 
reason is that the Newtonian part of the ranging delay (pQ) contains time in 
its left side, which must be transformed from one frame to another with tak- 
ing into account the post-Newtonian correction: T =>- t + v Q x/c 2 . This was 
done by Bertotti et al [6J who obtained that the ranging time delay in the 
heliocentric and barycentric frames must be related by the simple equation 



, „ s 2GM , 
At = (l+7-fc./3 ) — —In 



r 2 + n + r 12 



r 2 + r x - r 12 



(5) 



where r± 2 — \r 2 — fi|, r± — \r\\, r 2 = \r 2 \, (3 Q = v Q /c, and k is a unit vector 
along the light ray from the emitter to the observer. Formula ([5]) was derived 
previously in our work [8] (see also [9] for the case 7 = 1). 

Bertotti et al 0, section 4] noticed that our paper [S] neglected the post- 
Newtonian correction in the transformation of the time coordinate: T =^> 
t + Vqx/c 2 . Hence, [5] believed that [5] missed the velocity-dependent term 
in front of the logarithmic function in equation (jSJ). However, our paper [5] 
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dealt with the gravitational Doppler shift of the Cassini radio frequency, that 
is with the time derivative of the original heliocentric equation ([I]). Transfor- 
mation of the gravitational Doppler shift does not require to transform time in 
order to get all linear velocity-dependent corrections, because the heliocentric 
equation for the Doppler shift is already proportional to velocities of observer 
and emitter. Hence, the only transformation, which remains to complete is 
the transformation of the velocities, which can be done after making use of 
transformations ([3]) of spatial coordinates differentiated with respect to time 
(see [5J equations 12-15]). Our equation for the gravitational Doppler shift, 
z gr , written down in the bary centric coordinates is [5]I8"] 

V = _ ( — & + — • ol b - (1—P2 + I—P1 ~ As) • oc B , (6) 
Vri2 ri2 / V r 12 r 12 J 

where (3\ = (l/c)dxi/dt - velocity of the emitter, f3 2 = (l/c)dx 2 /dt - velocity 
of the observer, 7 = 7 — 1, and the dot between two vectors denotes a scalar 
dot product. Vector 

olb = a &~^ d ' ( 7 ) 

where R & - radius of the Sun, d - radius-vector of the impact parameter of 
the radio signal with respect to the Sun, d = \d\, and a© = 8.5 x 10 -6 rad is 
the solar gravitational deflection of light on its limb. 

The first term in the right side of equation ([6]) was obtained by Bertotti and 
Giampieri [10] under condition that the Sun does not move. The second term 
in the right side of equation ([6]) would lead to deviation from general relativity 
in case of 7 7^ 0. The solar velocity (3 & also enters this term, has the same 
magnitude as the terms with 7 and, hence, correlates with the measurement 
of 7. This olb(3q ~ 3 x 10~ 13 correction comes from the time derivative of 
the argument of the logarithmic function in equation (jHJ) as shown in [5]. 
Differentiation of the solar-velocity term in front of the logarithm in equation 
(jSJ) would give corrections of the next order of magnitude to the right side of 
the Doppler shit equation (jHJ). We did not forget that term but neglected it 
due to its smallness. We conclude that the remark Bertotti et al [6] on that we 
missed the velocity term in front of the logarithm in equation (jSJ) is irrelevant 
for our Doppler shift formula calculations [5]. 

The paper by Bertotti et al [6j also claims that the velocity-dependent terms 
appear in the time delay (JSj) only in front of the logarithmic function. This 
claim is rather naive and was the reason for misinterpretation of the impact 
of the solar motion on the Cassini data present in paper [B] . It only seems like 
that the argument of the logarithm in equation §5$) does not contain velocity- 
dependent terms explicitly. One should keep in mind that the distance r 12 is 
not taken on a single hypersurface of constant time but connects two different 
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positions of the Sun, which are not the same in the barycentric coordinates 
because the Sun moves. The post-Newtonian expansion of distance r 12 yields 



r\2 







owl) , 



(8) 



where r = \r\ and r = x 2 — X\ = kr is a null-cone vector connecting the 
observer to the emitter along the radio wave path from the Cassini to the 
observer. The post-Newtonian expansion of the ranging delay (j5J is 



„ ,2GM, 
At = (l+ 7 -fe./3 )^-lB 
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which explicitly reveals the presence of the velocity-dependent terms in the 
argument of the ranging time delay. Equation (Q has been derived in our paper 
[5J, and its logarithmic part is just a partial derivative of the heliocentric time 
delay with respect to the PPN parameter 7 as shown in [5J equation 25]. 

Bertotti et al. [B] claimed that expression for the ranging time delay is not 
used in the ODP code and can not be applied for theoretical analysis of the 
Cassini experiment as we did in [5]. Therefore, Bertotti et al. [6] have concluded 
that our numerical estimates of the gravitational shift of frequency caused by 
motion of the Sun with respect to the barycenter of the solar system as given in 
[5], are incorrect. These statements of Bertotti et al. [5] can not be accepted by 
any rationally thinking researcher as the authors of [6 J have trivially overlooked 
that expression (Q has exactly the same logarithmic function as in equation 
(J3J) with the argument expressed in terms of the null-cone distance r and 
velocity of the Sun, v Q , which are related to distance r\ 2 via self-consistent 
mathematical transformation (jSJ). The equivalence of equations (jSJ and 
means that our numerical estimates and theoretical conclusions given in paper 
[5J with regard to the impact of the solar motion on the Cassini measurement 
of 7 parameter, are firmly confirmed by the independent study of Bertotti et 
al 0. 
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